Linear response of a grafted semifiexible polymer to a uniform force field 
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We use the worm-like chain model to analytically calculate the linear response of a grafted semi- 
flexible polymer to a uniform force field. The result is a function of the bending stiffness, the 
temperature, the total contour length, and the orientation of the field with respect to that of the 
grafted end. We also study the linear response of a worm-like chain with a periodic alternating 
sequence of positive and negative charges. This can be considered as a model for a polyampholyte 
with intrinsic bending siffness and negligible intramolecular interactions. We show how the finite 
intrinsic persistence length affects the linear response to the external field. 
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I. INTRODUCTION 

Semifiexible polymers are macromolecules with a bend- 
ing stiffness intermediate between that of an absolutely 
flexible (Gaussian) chain and a rigid rod. A measure of 
the directedness of a polymer is its persistence length, 
which is proportional to the bending stiffness. Semifiex- 
ible polymers have a presistence length of the order of 
the total contour length. They have been the subject 
of extensive theoretical and experimental study in re- 
cent years primarily because many biologically important 
macromolecules fall in this class. The structural elements 
of the cytoskeleton (actin filaments, microtubules, inter- 
mediate filaments) and double-helix DNA exhibit elastic 
behavior dominated by their bending stiffness^. 

Single-molecule experiments have allowed the study 
of the influence of external forces or force fields on the 
conformational properties of semifiexible polymers 3-4-5-6 . 
Their theoretical study is considerably more challenging 
compared to that of flexible chains and one of the main 
reasons is the lack of scaling properties except for limiting 
cases (e.g., weakly bending polymers). In the worm- like 
chain model, pulling a semifiexible polymer at its ends 
is formally analogous to the Stark effect of a quantum 
rotator 7 . This analogy has led to semi-analytical solu- 
tions of the force-extension problem. The response to 
a force field is even more complicated as the quantum 
analogy involves a non-local (in imaginary time) inter- 
action. Further complications arise from intramolecular 
interactions between different polymer segments which 
may become important in realistic experimental situa- 
tions where a polymer is subject to a hydrodynamic flow 
or it is charged and placed in an electric field. 

Marko and Siggia- have derived an approximate ana- 
lytical field-extension relation for a semifiexible polymer 
which has one end fixed and is free to rotate about it. 
Their result for strong fields appears to be in good agree- 
ment with an experiment done with DNA in an external 
electric field^. Lamura et al*£> have calculated confor- 
mational properties of a grafted semifiexible polymer in 
a uniform force field in two dimensions using recursion 
relations in the weakly bending approximation which is 



valid at strong fields and/or large persistence lengths. 
We should mention that the standard Gaussian polymer 
model is inadequate to describe stretching in strong fields 
because it is infinitely extensibl e) 11 ' 12 . In contrast, the 
worm-like chain is characterized by a fixed total contour 
length constraint and this problem is avoided. 

Charged polymers with two types of charge (positive 
and negative) along their backbone are called polyam- 
pholytes. They have attracted a lot of attention because 
of their inherent theoretical interest and their relevance 
to the protein folding problem^ 3 -. There have been studies 
of their response to an electric fieldi^i 5 ^ 6 ^ 7 - but most of 
the previous works deal with their conformational prop- 
erties due to the intramolecular interactions in the ab- 
sence of an external field. They all consider charge dis- 
tributions on backbones described as flexible or freely- 
jointed chains and do not take into account the possibility 
of an internal bending stiffness. 

In the present paper we calculate the linear response 
of a grafted worm-like chain to a uniform electric field 
for arbitrary bending stiffness and field orientation. In 
Sec. II we describe the model. In Sec. Ill we calculate 
the linear response of the orientation and the extension 
of a uniformly charged worm-like chain and compare the 
results with the response to a force applied at the free 
end. fn Sec. IV we calculate the linear response of the 
extension of a periodic alternating polyampholyte with 
internal bending stiffness. 



II. THE MODEL 

The worm-like chain is a fluctuating, continuous, lo- 
cally inextensible line with a fixed total contour length 
L. In the absence of any external force, its effective free 
energy functional depends only on the bending (curva- 
ture) and is given by 



Wo[{r(s)}] = - 



ds 



dt(s) 



ds 



(2.1) 



where t(s) = dr(s)/ds is the unit tangent vector of the 
curve r(s) at arc length s, and k is the bending stiffness^. 
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The correlation length of the unit tangent vector along 
the polymer contour is the persistence length, L p , which 
is related to the bending stiffness via L p — 2K/[ksT(d — 
1)], where d is the dimensionality of the embedding space. 
Throughout this paper, we consider a worm-like chain 
which is grafted at s = 0, that is, both the position of this 
point and the orientation of the related tangent vector are 
kept fixed. 

The interaction with a uniform external field E is ex- 
pressed by adding to Ho the term 

Hj[{r( B )}] = -E. [ L dsX(s)[r(s)-v(0)} 
Jo 

= -E • f dsX(s) [ ds't(s') , (2.2) 
Jo Jo 

where A(s) is a phenomenological linear charge density. 
This simple model neglects the intramolecular interaction 
between different segments of the same polymer and also 
the steric self- avoidance. The latter is expected to be 
negligible for L p > L. For polymers with a large aspect 
ratio, like double-helix DNA or F-Actin, self-avoidance 
effects become important only for L 3> L p . 

If 0(s) is the angle between the tangent vector t(s) and 
the field E, we can write 



H = H -E dsX(s) ds'cosO(s') 



(2.3) 



where the first average in the rhs of the equation is taken 
over the Boltzmann weight associated with the energy 
given in Eq. 12.3fl and the second average is taken over 
the Boltzmann weight associated with the energy given 
in Eq. I|2.1|l . To lowest order in E, 



6(cos9(s)} 



EX 
k^f 

L 



ds' / ds" (cos 8(a) cos 0(s"))o 



-(cos0(s)) o / ds' / ds" (cos 9 (s"))o| • (3.2) 
Jo Jo ' 

In this Section, we assume that the linear charge den- 
sity A is constant along the polymer contour. Using the 
prescription of Eq. (|2.4H , one obtains 



(cos0(s))o = cos0(O) exp(— s/£ p ) 



(3.3) 



and 



(cos6(s) cos0(s'))o = cxp[— (s — s')/L p ] 

x { ~ + exp(-M/L p ) [ cos 2 0(0) - ~] } , (3.4) 

where s > s' . Substituting Eqs. I|3.3|) and l|3.4|l into Eq. 
(|3.2|) and breaking the domain of integration of the first 
double integral in the rhs of the equation into parts with 
a well-defined arc-length ordering, we obtain the final 
result: 



The orientational probability distribution of a worm- 
like chain having a conformation with a tangent vector 
t(0) = to at one end and a tangent vector t(L) = at 
the other is given by a path integral which is formally 
analogous to the density matrix element, in the angle 
representation, of a quantum rigid rotator—. We denote 
it by G(ti,i|to, 0). It can be calculated analytically and 
it has a spectral representation in terms of spherical har- 
monics. In 3 dimensions, we can integrate out the az- 
imuthal angle to obtain G[9(s), s\8(s'), s']. Thus we are 
able to calculate orientational correlations using 

(cos0(s„)...cos0(si)) o 

= / d9(L)sm9{L) [ d9(s n ) sin 9 (s n )... [ d9( Sl ) sm9{s x ) 



S(cos9(s)) 



EXL 
~~k^T 



2 s 2 L 2 L 



exp(-s/L p ) 



+i exp[(s - L)/L p ] - i exp[-(s + L)/L p ] 

+ (cos 2 0(O)-i) [Ii-exp(-3a/£p) 
3 

+ exp[-(2s + L)/L P ] - ~ exp(-3s/L p ) 
3 

+ 2 exp(-s/L p ) - exp[-(s + L)/L p ] 



~—exp(-s/L p )-~-?-exp(-3s/L p )]) . (3.5) 



xG[6(L),L\d(s n ),s n ] cosd(s n )G[d( 

^n)i Sri 

x...cos0(si)G[0(si),si|0(O),O] , 



where L > s n > s n _i > ... > si > 
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-i))S n _i] jg intstructive to consider the flexible and the 

(2.4) weakly bending limits of this result. In the flexible limit 
where L 3> L p , the anisotropy associated with the 0(0)- 
dependence of the response drops out and we get 



III. LINEAR RESPONSE OF A UNIFORMLY 
CHARGED SEMIFLEXIBLE POLYMER 

We define the change in the orientation of the polymer 
due to the applied force field by 



5{cos9(s)) = (cos9(s))e — (cos0(s))o 



(3.1) 



5(cos9(L)) 



3k R T 



+ o 



— ) exp ( - — 

lj p L, p 



(3.6) 



The leading term is precisely the response of a freely 
rotating dipole moment /i = AL 2 to a weak field 
E . The average polarization of such a system is 
(/zcos0) = (l/2)/d0sin0//cos0exp[-/i£cos0/(fc s T)]. 
On the other hand, in the weakly bending limit where 
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polymer contour: 




1e2 .1e3 .1e4 



FIG. 1: Rescaled linear response coefficient for the extension, 
C" = 5 i?||fcsT / (EXLpL 2 ) , as a function of the chain stiffness, 
I = L/Lp, for 8q — it/ 2 (upper curve) and do — (lower 
curve) . 



L <c L p , the response strongly depends on 6(0). For 
0(0)^(0,70, 



<5(cos0(L)) = ^sin 2 0(O) + o((A) 4 ) . (3.7) 



For a purely bending field (0(0) = tt/2), we recover the 
mechanical bending of a rod with siffness k which is 
known from the classical theory of elasticity^. For an 
elongational (or compressional) field with 8(0) = (0, 7r), 



6 (cos 6 (L)) 



EXL 4 k B T 
12k? 



+ ((i-)=). (3 . 8) 



We note that in this case the linear response becomes 
proportional to the temperature and it vanishes at zero 
temperature. This vanishing is an indication of the ex- 
istence of a finite critical field above which the compres- 
sional deformation becomes unstable as predicted by the 
classical theory of elasticity of rods (buckling to a field 

An experimentally more accessible quantity is the re- 
sponse of the end-to-end vector in the direction of the 
field: 



5R l{ = ((R) E - (R) ) 



E 

lEi 



(3.9) 



where 



R = f dst(s) . (3.10) 
Jo 

We easily obtain it by integrating Eq. I|3.5|l over the 



-— exp(-L/L„)] + ~ [(A) 2 + cxp(-2L/L p ) 

Lip 6 Lip 



+36] 
5 



1, 



'0(0)--)[ 



1 L 
3 L r . 



4 7 

4 exp(-L/L p ) + - - — exp(~3L/L p ) 
cxp(-2L/L p ) + i A cxp(-L/L p )] } . (3.11) 



As with the response of the orientation, we gain insight 
by considering the limiting cases of flexible and weakly 
bending polymers. In the flexible limit (L 3> L p ), 



SR11 = EXL ^ 



(3.12) 



The leading term is the result that one obtains for the 
response to a uniform field of a Gaussian chain with Kuhn 
length equal to 2L p . Obviously, this result cannot be 
valid in the strong field regime because the ~ I? scaling 
would become incompatible with the local inextensibility 
of the polymer. Indeed, the response of a worm-like or a 
freely-jointed chain to a strong field has different scaling 
behavior and is free from this pathology as it is shown 
in Refs^fii, andi^. For weak fields, however, there is no 
inconsistency between Eq. H3.12|l and the constraint of 
fixed contour length. In the weakly bending limit (L <c 
L p ), 



SRn = 



EXL 4 



■ sin 



0(O) + O((A) 5 ) (3.13) 



for 9(0) ^ (0,tt), and 



SRn 



EXk B TL 5 
20^2 



(3.14) 



for 6(0) = (0, tt). As before, the bending response in the 
weakly bending limit reproduces the result from the clas- 
sical elasticity of rodsiS and the temperature dependence 
in Eq. 13.14fl is a sign of its entropic origin. The crossover 
between the two scaling limits as the chain stiffness varies 
is illustrated in Fig. 1. 

In order to compare the linear response to a field to 
the linear response to a force exerted at the free end, it 
is instructive to define an "effective point charge." It is a 
point charge placed at the free end which, subject to the 
field E, would yield the same response as that due to the 
total charge XL which is uniformly distributed along the 
polymer contour. In Fig. 2, we plot the "effective point 
charge" (in units of XL) as a function of L/L p for the re- 
sponse of the extension, SRn , (upper curve) and for that 
of the orientation, S (cos 6(L)), (lower curve). The equa- 
tions giving the response to a force applied at the free 
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FIG. 2: Effective point charge for the response of the ex- 
tension (upper curve) and for the response of the orientation 
(lower curve) as a function of I = L/L p for 8q — 0. 




FIG. 3: Linear response coefficient for the extension of a pe- 
riodic polyampholyte, C = 5R\\kBT/(EXoLp), as a function 
of I = L/Lp, for 0o = vr/2 and l m /L p = 1/20. 



end of a grafted worm-like chain are cited in Appendix 
A. We note that the effective charge for the extension 
tends to (1/2) XL in the flexible limit which is just the 
average charge along the polymer contour. In the weakly 
bending limit, it decreases to (3/10)AL for 0(0) — (0,ir) 
and to (3/8) XL for any other orientation. 



IV. ALTERNATING POLYAMPHOLYTES 

In this Section, we extend our study of a charged worm- 
like chain in an external electric field to a grafted semi- 
flexible polymer with a periodic alternating charge se- 
quence along its contour. We model such a sequence by 
a sinusoidal linear charge density: 



X(s) = Aq sin(fcs) , 



(4.1) 



where k = 2ir/l m with l m being the arc length of the 
"elementary dipole." This model is analytically tractable 
and the linear response of the extension formally reads: 

SR ]{ = 

— [ ds{ [ ds'sm(ks') [ ds" (cos 0(a) cos 6(s")} 
k B T J I J Jo 

-(cos6»(s)} ^ ds'sin(fcs') J ds"(cos 6>(s"))o} . (4.2) 

The explicit final result appears too long and cumber- 
some and we report it in Appendix B. Here, we only 
highlight its main features. The sign of the response os- 
cillates depending on the sign of the charge at the tip of 
the chain (Fig. 3). Specifically, SRii will have the same 
sign as that of the charge in the last polymer segment 
of arc length i m /2. This is a manifestation of the same 
phenomenon as the "odd-even effect" discussed in Refs^ 



andi& for free (non-grafted) alternating polyampholytes. 
In those papers it is shown that a flexible polyampholyte 
in an external field stretches when the total number of 
charges is odd, while it collapses when the total number 
of charges is even (zero net charge) . In the case consid- 
ered here, because of the different boundary conditions, 
the polymer deforms in the direction of the external field 
when the charge at the last segment of length l m /2 is pos- 
itive and it deforms in the opposite direction when this 
charge is negative. The remarkable point is the strong 
dependence of the response on the charge of the last seg- 
ment despite the very large number of segments. 

In order to understand the interplay between the three 
characteristic lengths, L,L pi and l m , we focus on the 
response of a chain with total contour length L — ml m , 
where m (the number of dipoles) is an integer. In this 
case, as well as when L = (m + l/2)l m , the response is 
peaked. The result reads: 

EXoLl 

R]l ~ k B T6(9 + K 2 )(l + K 2 )(A + K 2 )K 2 
x \qK 7 - 87rif 6 TO - 288ttto + 2K 7 exp(~2L/L p ) 

-8K 7 cxp(-L/L p ) + 72K 3 exp(-2L/L p ) 
-104iv 5 cxp(-i/ J L p ) - 288i i 3 exp(~ J L/L p ) 
+216K 3 + 26K 5 exp(-2L/L p ) + 7SK 5 - 3927r!T 2 m 

-1127rA' 4 m + ( cos 2 0(0) - -) [78A 5 exp(~2L/L p ) 

3 

-297A 3 exp(-L/Lp) - 23A 3 exp(-3L/L p ) 
-87 'A 5 exp(- L/Lp) - 25A 5 exp(-3L/L p ) 
+216A 3 cxp(-2L/ J L p ) + 6A 7 exp(-2L/Lp) 
-2A 7 exp(-3 J L/L p ) - 6K 7 exp(-L/L p ) 

+2K 7 + 34A 5 + 104A 3 ] j , (4.3) 
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where K = 2n L p /l m . 

In the flexible limit (L ^> L p ) 



SR \\= , T 

k R T 



ZK 2 



we obtain: 



K 



3(1 + K 2 )(9 + K 2 ) 



[27 + 3K 2 + ( cos 2 0(0) - -) (13 + K 2 ) 
L 



(4.4) 



The leading term is the response of a neutral Gaussian 
chain with a point charge at its free end equal to the 
average charge in the last segment of arc length l m /2 of 
the polyampholyte, which is — (l/2)Aol m /ir. The correc- 
tion from the charge in the bulk depends on the ratio 
l m /L p and it is maximal when l m / ' L p w 2tt. We have to 
distinguish between the case where the chain is flexible 
on the scale of the elementary dipole, that is, L p -C l m , 
and the case where the chain is weakly bending on that 
scale. To leading order in L/L p , the relative difference 
between elongational and bending response defined as 
[5R {1 (9 = 0)-5R ll (0 = n/2)]/6R\\(9 = 0) is (2/3)L p /L 
for L p > l m , and ~ (L p /l m ) 2 (L p /L) for L p <C l m . We 
point out the qualitative difference between the two cases 
which could be used as a way to experimentally probe the 
polyampholyte stiffnes on the scale of l m . 

In the weakly bending limit (L <C L p ), for 9(0) ^ 

(0,7T), 



SRu 



O 



^^{Vm 3 + 2.m}sin^(0) + 



-)" 



whereas for 9(0) = (0, 7r), 
_ EX Ll 8 



4_4 

— 7T m 



O 



(4.5) 



(4.6) 



For large m, as in the flexible limit, we recover the re- 
sponse of a neutral worm-like chain with a point charge 
at its free end equal to the average charge in the last 
segment of arc length l m /2 of the polyampholyte. Notice 
that in the case of elongational or compressional fields, 
we do not need to assume large m in order to get this 
cancellation of the contribution from the charges in the 
bulk of the chain. 

If we modify the linear charge density of Eq. 14. 1|) by 
introducing a phase, that is, A(s) = Ao sin(fcs + </>), the 
results presented in this Section will remain unchanged 
up to a phase shift. For example, if we have <j) = 7r/2, the 
response will be peaked when the total contour length is 
L = (m + l/4)/ m or (m + 3/4)i m , where m is an integer. 



V. CONCLUSIONS 

In this paper we have calculated analytically the lin- 
ear response of a grafted worm-like chain to a uniform 
force field. The response assumes scaling forms in the 
limiting cases of weakly bending and flexible polymers 
and we have obtained explicit results for polymers of ar- 
bitrary stiffness. We have discussed how the response 
to a field differs from the response to a force exerted at 
the free end. We have considered a uniformly charged 
chain and a periodic alternating polyampholyte and we 
have demonstrated the strong dependence of the response 
on the distribution of charge along the polymer contour. 
In the latter case, we have shown how the interplay be- 
tween the internal persistence length and the characteris- 
tic length of the charge modulation affects the response. 



APPENDIX A 

For the sake of completeness and in order to facili- 
tate comparisons, in this Appendix, we cite the main 
results concerning the linear response of a grafted worm- 
like chain to a force applied at its free end 2 ^. In place of 
Eq. 1)2.2(1 . the interaction with a force F is expressed by 
adding to TL a the term 



TLf = — F 



dst(s) . 



(Al) 



Using the orientational correlations of the free chain, one 
obtains the linear response of the tip orientation, 

F L ( 1 2 1 
5<cos0(L)} = T -*U - xexp(-L/L p ) + - exp(-2L/L p ) 



k B Tl3 3 
-(cos 2 9(0) - i)[exp(-2L/L p ) 

\ exp(-L/L p ) - ~ exp(-3L/L p )] j 



and the linear response of the extension, 



(A2) 



FLz 



2L 4 

-1+ 3^ + ^ C M-L/L p ) 



~ exp(-2L/L p ) + ( cos 2 9(0) - 1) [ - i 
1 



exp(— 3L/L p ) + exp(—L/L p ) 
o 

- exp(-2L/L P ) - exp(-2L/L p )] } 



(A3) 
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APPENDIX B 

In this Appendix, we report the explicit result for the linear response of a periodic alternating polyampholyte of 
arbitrary length as obtained from Eq. i|4.2[l . 

6R *=- k B T6(9 + K^1°S)(4 + K*)K* eXP ^) cos(AZ) + ^ exp(-Z) cos(AZ) 

+340A 3 exp(-Z) cos(AZ) - 144 sin(AZ) - 2A 6 exp(-2Z) sin(AZ) - 26A 4 exp(-2Z) sin(AZ) 
-2K 7 exp(-2Z) cos(Kl) + 4K 5 - 52 A 3 exp(-Z) - 30 sin( AZ) A 4 - 144A exp(-Z) - 4A 5 exp(-Z) 
+144A + 4 cos(Kl)K 7 l + 56 cos( AZ)/fc 5 Z + 144 cos(Kl)Kl - 26A 5 exp(-2Z) cos( AZ) + 196 cos(AZ)A 3 Z 
+144 exp(-Z) cos(AZ) A - 124A 2 sin(Kl) ~ 2A 6 sin(Kl) - 6 cos(AZ) A 7 - 268 cos(Kl)K 3 - 82 cos(Kl)K 5 
+52K 3 + 52 exp(-0 sin( ATZ) X 4 - 144 cos(Kl)K + 144 exp(-Z) sin(A'/) K 2 ~ 72 A 3 exp(-2Z) cos(Kl) 

-72K 2 exp(-2l) sin(Kl) + 4A 6 exp(-J) sin(A0 + ( cos 2 6(0) - \) [81 A 5 exp(-l) cos(AZ) 

+255 exp(-Z) A 3 cos(AZ) + 6 exp(-Z) A 7 cos(AZ) + 23 exp(-3Z) A 3 cos(AZ) + 108 cos(AZ) cxp(-Z) A 
+25 exp(-3Z) cos(AZ) A 5 + 45 exp(-3Z) A 4 sin(AZ) + 2 exp(-3Z) A 7 cos(AZ) + 3 exp(-3Z) A 6 sin(AZ) 
+42 exp(-3Z) A 2 sin(AZ) - 216A 3 exp(-2Z) - 78A 5 exp(-2Z) cos(AZ) - 6A 6 exp(-2Z) sin(AZ) 
-6A 7 exp(-2Z) cos( AZ) + 6 exp(-Z) A 5 + 42 exp(-Z) A 3 - 108 exp(-Z) A + 39 A 4 exp(-Z) sin(AZ) 
-98 cos(AZ) A 3 - 28 cos(AZ) A 5 + 108 exp(-Z) A 2 sin(AZ) + 3 exp(-Z) A 6 sin(AZ) - 72 cos(AZ) A 
+72A - 6A 5 + 72A - 6A 5 + 6A 3 - 2 cos(AZ) A 7 - 216A 2 exp(-2Z) sin(AZ) 

-78A 4 cxp(-2Z)sin(AZ)]| , (Bl) 
where Z = L/L p . 
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